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Abstract 

In geophysics and seismology, it is a common knowledge that the quality factors Q of 
the mantle and crust materials scale as the tidal frequency to a positive fractional power 
(Karato 2007, Efroimsky & Lainey 2007). In astronomy, there exists an equally common 
belief that such rheological models introduce discontinuities into the equations and thus 
are unrealistic at low frequencies. We demonstrate that, while such models indeed make 
the conventional expressions for the tidal torque diverge for vanishing frequencies, the 
emerging infinities reveal not the impossible nature of one or another rheology, but a 
subtle flaw in the underlying mathematical model of friction. Flawed is the common 
misassumption that damping merely provides phase lags to the terms of the Fourier 
series for the tidal potential. A careful hydrodynamical treatment by Sir George Darwin 
(1879), with viscosity explicitly included, had demonstrated that the magnitudes of the 
terms, too, get changed - a fine detail later neglected as "irrelevant". Reinstating of 
this detail tames the fake infinities and rehabilitates the "impossible" scaling law (which 
happens to be the actual law the terrestrial planets obey at low frequencies). 

This preprint is a pilot paper. A more comprehensive treatise on tidal torques is to 
be published. (Efroimsky & Williams 2009) 
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1 Introduction. 



We are considering bodily tides in a primary perturbed by a secondary. Each elementary 
volume of the primary is subject to a tide-raising potential, which in general is not periodic 
but can be expanded into a sum of periodic terms. We shall assume that the primary is 
homogeneous and incompressible. Although simplistic, the model provides a good qualitative 
understanding of tidal evolution of both the primary's spin and the secondary's orbit 



2 Linearity of the tide 
2.1 Two aspects of linearity 

We assert deformations to be linear. Each tidal harmonic Wi of the potential disturbance 
produced by the secondary generates a linear deformation of the primary's shape, while each 
such deformation amends the potential of the primary with an addition proportional to the 

Love number ki . Linearity also implies that the energy attenuation rate {E (x)) each 
frequency x depends solely on the value of x and on the amplitude Epeakix) ? and is not 
influenced by the other harmonics. Thus, 

{Eix)) = -^^§^ (1) 

or, equivalently: 

^E^ycUx) = - ^^^1^ , (2) 

AEcycieix) being the one-cycle energy loss. The so-defined quality factor Q{x) corresponding 
to some frequency x is interconnected with the phase lag e{x) corresponding to the same 
frequency Q 

Q = cot |e| . (3) 

This makes the essence of the standard, linear theory of bodily tides. The model permits for the 
freedom of choice of the functional dependency of the quality factor upon the tidal frequency. 
Whatever the form of this dependence, the basic idea of the theory is the following: the tide- 
raising potential is expanded over periodic terms, whereafter (a) the material's response is 
assumed to be linear, and (b) the overall attenuation rate is asserted to be a sum of rates 
corresponding to the involved frequencies. Leaving the first of these assertions intact, we shall 
make in the subsequent subsection an aside on the applicability limits of the second assertion. 
In the Appendix below, we shall also probe some of the applicability limits of ([3]). 



2.2 Goldreich's admonition: 

a general difficulty stemming from nonlinearity 

Introduced empirically as a means to figleaf our lack of knowledge of the attenuation process in 
its full complexity, the notion of Q has proven to be practical due to its smooth and universal 

^ The interconnection becomes more complicated than Q when the values of Q are about ten or less - see 
the Appendix below. 
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dependence upon the frequency and temperature. At the same time, this empirical treatment 
has its predicaments and hmitations. Its major inborn defect was brought to hght by Peter 
Goldreich who pointed out that the attenuation rate at a particular frequency depends not only 
upon the appropriate Fourier component of the stress, but also upon the overall stress. This 
happens because for real minerals each quality factor Q{Xi) bears dependence not only on the 
frequency Xi itself, but also on the magnitude of the Xi~component of the stress and, most 
importantly, also on the overall stress. This, often-neglected, manifestation of nonlinearity 
may be tolerated only when the amplitudes of different harmonics of stress are comparable. 
However, when the amplitude of the principal mode is orders of magnitude higher than that of 
the harmonics (tides being the case), then the principal mode will, through this nonlinearity, 
make questionable our entire ability to decompose the overall attenuation into a sum over 
frequencies. Stated differently, the quality factors corresponding to the weak harmonics will 
no longer be well defined physical parameters. 

Here follows a quotation from Goldreich (1963): 
"... Darwin and Jeffreys both wrote the tide-raising potential as the sum of periodic potentials. 
They then proceeded to consider the response of the planet to each of the potentials separately. 
At first glance this might seem proper since the tidal strains are very small and should add 
linearly. The stumbling block in this procedure, however, is the amplitude dependence of the 
specific dissipation function. In the case of the Earth, it has been shown by direct measurement 
that Q varies by an order of magnitude if we compare the tide of frequency 2uj — 2n with 

o 

the tides of frequencies 2uj — n , 2uj ~ 3n , and . This is because these latter tides have 
amplitudes which are smaller than the principle tide ( of frequency 2uj — 2n ) by a factor of 
eccentricity or about 0. 05. It may still appear that we can allow for this amplitude dependence 
of Q merely by adopting an amplitude dependence for the phase lags of the different tides. 
Unfortunately, this is really not sufficient since a tide of small amplitude will have a phase lag 
which increases when its peak is reinforcing the peak of the tide of the major amplitude. This 
non-linear behaviour cannot be treated in detail since very little is known about the response of 
the planets to tidal forces, except for the Earth. " 

On these grounds, Goldreich concluded the paragraph with an important warning that we 
"use the language of linear tidal theory, but we must keep in mind that our numbers are really 
only parametric fits to a non-linear problem. " 

In order to mark the line beyond which this caveat cannot be ignored, let us first of all 
recall that the linear approximation remains applicable insofar as the strains do not approach 
the nonlinearity threshold, which for most minerals is of order 10~® . On approach to that 
threshold, the quality factors may become dependent upon the strain magnitude. In other 
words, in an attempt to extend the expansion (P - [2]) to the nonlinear case, we shall have 
to introduce, instead of Qixd , some new functions Q{Xi, -£^peafc(xj> ) • (Another 

complication is that in the nonlinear regime new frequencies will be generated, but we shall 
not go there.) Now consider a superposition of two forcing stresses - one at the frequency Xi 
and another at X2 ■ Let the amplitude E {xi ) be close or above the nonlinearity threshold, 
and E ^^^{x2) be by an order or two of magnitude smaller than E ^^^(Xi) • To adapt the 
linear machinery (P - [2]) to the nonlinear situation, we have to write it as 

(E) = {E,} + {E2) = - X^ o(f'7^'''].)) - ^2 nrv F \v^f V )) ' 

the second quality factor bearing a dependence not only upon the frequency X2 and the 
appropriate magnitude E^^^^{x2) ' tiut also upon the magnitude of the first mode, E^^^i_{xi) 5 
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- this happens because it is the first mode which makes a leading contribution into the overall 
stress. Even if (jl]) can be validated as an extension of ([1] - [2]) to nonlinear regimes, we should 
remember that the second term in (jlj) is much smaller than the first one (because we agreed 
that E^^^^{x2) ^ ^peakiXi) )• This results in two quandaries. The first one (not mentioned 
by Goldreich) is that a nonlinearity-caused non-smooth behaviour of Q{Xi ^ E ^^f^{Xi) ) will 
cause variations of the first term in (jl]), which may exceed or be comparable to the entire 
second term. The second one (mentioned in the afore quoted passage from Goldreich) is the 
phenomenon of nonlinear superposition, i.e., the fact that the smaller-amplitude tidal harmonic 
has a higher dissipation rate (and, therefore, a larger phase lag) whenever the peak of this 
harmonic is reinforcing the peak of the principal mode. Under all these circumstances, fitting 
experimental data to (jlj) will become a risky business. Specifically, it will become impossible 
to reliably measure the frequency dependence of the second quality factor; therefore the entire 
notion of the quality factor will, in regard to the second frequency, become badly defined. 

The admonition by Goldreich had been ignored until an alarm sounded. This happened 
when the Jet Propulsion Laboratory's Lunar-ranging team applied the linear approach to 
determining the frequency-dependence of the Lunar quality factor, Q{x) ■ They obtained (or, 
as Peter Goldreich rightly said, fitted their data to) the dependency Q{x) ~ (Williams et 
al. 2001). The value of the exponential for the Moon turned out to be negative: a = — 0.07 , 
a result firmly tabooed by the condensed-matter physics for this range of frequencies (Karato 
2007). 

One possible approach to explaining this result may be the following. Let us begin with a 
very crude estimate for the tidal strains. For the Moon, the tidal displacements are of order 
0.1m (zero to peak). The most rough estimate for the strain can be obtained through dividing 
the displacement by the radius of the body. While this ratio is almost twenty times less than 
the 10~^ nonlinearity threshold, we should keep in mind that in reality the distribution of the 
tidal strain is a steep function of the radius, with the strain getting its maximum near the 
centre of the body, as can be seen from equations (48.17) in Sokolnikoff (1956). The values 
of strain can vary in magnitude, over the radius, by about a factor of fivej^ Another order 
of magnitude, at least, will come from the fact that the Moon is inhomogeneous and that its 
warmer layers are far more elastic than its rigid surface. As a result, deep in the Lunar interior 
the tidal strain will exceed the afore mentioned nonlinearity threshold. Hence the uncertainties 
in determination of Q{x) by the Lunar-ranging team. 

We shall not dwell on this topic in quantitative detail, leaving it for a future work. Our 
only goal here has been to draw the readers' attention to the existing difficulty stemming from 
the shortcomings of the extension of (P - [2]) to nonlinear regimes. 

3 Darwin (1879) and Kaula (1964) 

The potential produced at point R = {R, X , (p) by a mass M* located at r* = (r*. A* , 0*) is 

oo ' I 

W{R,n = - ^^f^f^) ^|l^!^(2-5ojP,^(sin0)P,^(sinr) cosm(A-A*) . (5) 

^ 1=2 m=0 ^ +m,j. 

^ Accordingly, the tidal-energy radial distribution, too, is strongly inhomogeneous, with a maximum in the 
planet's centre - see Fig. 3 in Peale & Cassen (1978). 
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When a tide-raising secondary located at r * distorts the shape of the primary, the potential 
generated by the primary at some exterior point r gets changed. In the linear approximation, 
its variation is: 

1=2 V ^ / 

ki being the Ith. Love number, R now being the mean equatorial radius of the primary, 
R = {R, (p , X) being a surface point, r* = (r* , 0* , A*) being the coordinates of the tide- 
raising secondary, r = (r , , A) being an exterior point located above the surface point R 
at a radius r > R, and the longitudes being reckoned from a fixed meridian on the primary. 
Substitution of ([5]) into entails 

/?^'+^ ' (J- V 

U{r) = -GM*Y,ki Yl L ; (2-<^ojP/,n(sin0)P^^(sin0*) cosm(A-A*) . (7) 

1=2 r r* ^=oy^ + ^)- 

A different expression for the tidal potential was offered by Kaula (1961, 1964), who developed 
a powerful technique that enabled him to switch from the spherical coordinates to the Kepler 
elements (a*, e*, i*, Q*, uj*, JH* ) and (a, e, i, Q, u, A4) of the secondaries located at r * 
and r. Application of this technique to ([7]) results in 

Z=2 ^ ^ ^ ^ m=0 ^ ^ 

(8) 

Z oo i oo 

-Som) ^^Fimpii*) E Gipg{e*)'YFimh{i) 5Z ^'^i(^) COS -m6'*) - -m6')] 

p=0 q=—oo h=0 j=—oo 

where 

= {l-2p)uj* + {l-2p + q)M* + mQ* , (9) 

vimhj = {I -2h)uj + {I -2h + j)M + mn , (10) 

and 6' = 6'* is the sidereal angle. 

While ([8]) and ([7]) are equivalent for an idealised elastic planet with an instant response of 
the shape, the situation becomes more involved when dissipation-caused delays come into play. 
Kaula's expression ([H]), as well as its truncated, Darwin's version jl is capable of accommodating 
separate phase lags for each harmonic involved: 

1=2 ^ ^ ^ m=0 ^ ' 

(11) 

Z oo Z oo 

5om ) y^^Fimpji*) Glpq{e*)y^^Fimh{i) ^ Glhj{e) COS [{v*^pg - m9*) - {vimhj - m9) - eimpg] 

p=0 q=—oo h=0 j=—oo 

^ While the treatment by Kaula (1964) entails the infinite Fourier series ([5|), the developments by Darwin 
(1879) furnish its partial sum with \q\, \j\ < 2. For a simple introduction into Darwin's method see 
Ferraz-Mello, Rodriguez & Hussmann (2008). 
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where 



'Impq 



{I - 2p) u* + {l~2p + q)M* + m {Cl* - 9*) 



At 



Impq 



^Impq ^tlmpq 



± X*lmpq 



impq — -"Impq 



is the phase lag interconnected with the quahty factor via Qimpq = 
harmonic uj^^j^pq introduced in (fT2!) is 

cu*^^^ ^ {l-2p)dj* + {l-2p + q)M* + m{tl* - 0* 
while the positively-defined quantity 

I (/ - 2p) Cj* + {l-2p + q)M* + m {tl* 



cot e 



Impq I 



. The tidal 



(13) 



Xlmpq 



Impq 



9* 



(14) 



is the actual physical Impq frequency excited by the tide in the primary. The corresponding 
positively-defined time delay Atimpq depends on this physical frequency, the functional forms 
of this dependence being different for different materials. 

Formulae ([H]) and (fTTl) constitute the principal result of Kaula's theory of tides. Most impor- 
tantly, Kaula's formalism imposes no a priori constraint on the form of frequency- dependence 
of the lags. 



4 The Darwin-Kaula-Goldreich expansion 
for the tidal torque 

Now we are prepared to calculate the planet-perturbing tidal torque. Since in what follows we 
shall dwell on the low-inclination case, it will be sufficient to derive the torque's component 
orthogonal to the planetary equator: 



r 



- M 



dUjr] 
36 



(15) 



M being the mass of the tide-disturbed satellite, and the "minus" sign emerging due to our 
choice not of the astronomical but of the physical sign convention. Adoption of the latter 
convention implies the emergence of a "minus" sign in the expression for the potential of a 
point mass: —Gm/r. This "minus" sign then shows up on the right-hand sides of ([7]), 
([S]), and (ITTl) . It is then compensated by the "minus" sign standing in (ITH) . 

The right way of calculating dU (r ) / 89 is to take the derivative of ffTTl) with respect to 9 , 
and theiQ to get rid of the sidereal angle completely, by imposing the constraint 9* = 9 . This 
will yield: 

1=2 ^ ^ ^ ^ m=0 ^ ' p=0 q=-oo 



/i=0 



]=-ca 



Be mindful that our intention here is to differentiate not cos 



Impq 



-mO' 



(16) 



[vimh] - mB) 



but 



(^, 



■^{delayed) 
Impq 



me* 



{delayed) 



) ~ [vimh] - 



Hence tlie said sequence of operations. 
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In the case of the tide-raising satelhte coinciding with the tide-perturbed one, all the elements 
become identical to their counterparts with an asterisk. For a primary body not in a tidal lock 
with its satellitejfl it is sufficient to limit our consideration to the constant part of the torquej^ 
a part for which the indices (p, q) coincide with {h , j) , and therefore vimhj cancels with 
'"impq ■ This will give us: 

, ^ 2 2^ + 1 -2/-2 , ^ (/ — my. \ ^ 2 , ^ 2 

r = 2_^2 ki G M R a (I + rn)\ ^ ^ ^irnpi^) 2^ Gipgi^) sineimpq ■ (17) 

1=2 m=0 ^ p=0 q=-oo 

The expression gets considerably simplified if we restrict ourselves to the case of / = 2 . Since 
< m < I , and since m enters the expansion as a multiplier, we see that only m = 1 , 2 
actually matter. As < p < / , we are left with only six relevant F's, those corresponding 
to (Imp) = (210), (211), (212), (220), (221), and (222). By a direct inspection of the table of 
Fimp we find that five of these six functions happen to be 0{i) or 0( ) , the sixth one being 

^ ( 1 + cos i) = 3 + O(z^) . Thus we obtain, in the leading order of i : 



-t^220 

r = - V GM^ Ra-^G^ (e) /ta sine + O(zVQ) • (18) 

q= — oo 

The leading term of the expansion is 

r^,,, = \ GM'hR\-'sme,,,, . (19) 
Switching from the lags to quality factors via formulc0 

Qlmpq I cot €lfYipq \ , (20) 

we obtain: 

sgntu, sgncu, sgncj, 

sine, = sin|e, | sgncu, = ''""^ = = , ^"'"^ = -—^^ + OiQ-') , (21) 

Impq Impq Impq H. ^ 2 ^ / 2 O 

^'^ + ^Impq \ I + Q ^""Pl 

y Impq 



^ With a denoting the librating angle, the locking condition reads: 9 = fl + uj + Ai + 180° + a + 0{i^) . 
Insertion thereof into (fTS]) results in: i^;*„jpq = {I — 2p — m) uj* + {I — 2p + q — m) M* , where we have 
neglected —ma on account of a being extremely small. Clearly, the indices can assume more than one set of 
values corresponding to one tidal frequency. This way, the case of libration is more involved than that of tidal 
despinning, and deserves a separate consideration. 

^ The tide-raising and tidally-perturbed satellites being the same body does not yet mean that the indices 
{p , q) coincide with with {h , j) . These are two independent sets of indices, wherewith the terms of two Fourier 
series are numbered, expression (|16p being a product of those two series. This product contains a constant part, 
as well as short-period terms proportional to A4 and long-period terms proportional to uj . The short-period 
terms get averaged out over a period of the tidal flexure, while the long-period terms get averaged out over 
longer times, provided the periapse is precessing and not librating. Expression (jl7p furnishes the constant part 
of the torque. Fortunately, this is sufficient for our further calculations. 

^ The phase lag eimpq is introduced in (|12p . while the tidal harmonic cuimpq is given by ()13|) . The quality 
factor Qlmpq = \ cot eimpq I is, for physical reasons, positively defined. Hence the multiplier sgntoimpq in (|2ip . 
(As ever, the function sgn(x) is defined to assume the values -1-1 , —1 , or for positive, negative, or vanishing 
X , correspondingly.) 

Mind that no factor of two appears in (|20l - [2T|) , because e is a phase lag, not a geometric angle. 
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whence 



1=2 



- J2 GM' R a^' gI (e) k. 



Q 



q=-oo ■'220q 

Now, let us simplify the sign multiplier. If in expression f|T3|) for cu^^p^ we get rid of the 

redundant asterisks, replac^ M. with M.o + n ~ n, and set / = m = 2 andp = 0, the outcome 
will be: 

sgn uj^^^q = sgn 2ci;+(2 + g)n + 2fi— 26* = sgn 



c:; + 1 + 



n + n - e 



As the node and periapse precessions are slow, the above expression may be simplified to 



sgn 



1 + ^ 
2 



n 



9 



All in all, the approximation for the torque assumes the form: 

1 



r 



1=2 



- V G R a-^Gl^Je) Q ^ sgn 



- ] n — 



OiiyQ) + OiQ- 



That the sign of the right-hand side in the above formula is correct can be checked through the 
following obvious observation: for a sufficiently high spin rate 6 of the planet, the multiplier 

sgn (^^ + § ) n — 9 becomes negative. Thereby the overall expression for r acquires a 
"minus" sign, so that the torque points out in the direction of rotation opposite to the direction 
of increase of the sidereal angle 9 . This is exactly how it should be, because for a fixed q and 
a sufficiently fast spin the g's component of the tidal torque must be decelerating and driving 
the planet to synchronous rotation. 

Expansion (1221) was written down for the first time, without proof, by Goldreich & Peale 
(1966). A schematic proof was later offered by Dobrovolskis (2007). 



5 Can the quality factor scale as a positive power of the 
tidal frequency? 

As of now, the functional form of the dependence Q{x) for Jovian planets remains unknown. 
For terrestrial planets, the model Q ~ 1/x is definitely incompatible with the geophysical 
data. A convincing volume of measurements firmly witnesses that Q of the mantle scales as 
the tidal frequency to a positive fractional power: 

Q = , where a = 0.3 ± 0.1 , (23) 

® While in the undisturbed two-body setting A4 = Mo + n {t ~ to) and M = n , under perturbation these 
relations get altered. One possibility is to introduce (following Tisserand 1893) an osculating mean motion 
n{t) = ^/ a{t)^ , and to stick to this definition under perturbation. Then the mean anomaly will evolve as 
M = Mo{t) + n{t) dt , whence M = Mo{t) + n{t) . 

Other possibilities include introducing an apparent mean motion, i.e., defining n either as the mean-anomaly 
rate dA4 /dt , or as the mean- longitude rate dL/dt — dft /dt + duj/dt + dA4 /dt (as was done by Williams et 
al. 2001). While the first-order perturbations in a{t) and in the osculating mean motion -y/ fi/a{t)'^ do not have 
secular rates, the epoch terms typically do have secular rates. Hence the difference between the apparent mean 
motion defined as dL/dt (or as dA4/dt) and the osculating mean motion y/Ji/a(T)^. I am thankful to James 
G. Williams for drawing my attention to this circumstance (J.G. Williams, private communication). 
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8 being an integral rheological parameter with dimensions of time. This rheology is incompat- 
ible with the postulate of frequency-independent time-delay. Therefore an honest calculation 
should be based on averaging the Darwin-Kaula-Goldreich formula (1221) . with the actual scal- 
ing law (12511 inserted therein, and with the appropriate dependence At^^^^ ( Ximpq ) taken into 
account n 



5.1 The "paradox" 

Although among geophysicists the scaling law ( 123|) has long become common knowledge, in 
the astronomical community it is often met with prejudice. The prejudice stems from the fact 
that, in the expression for the torque, Q stands in the denominator: 

. ~ i . (24) 

At the instant of crossing the synchronous orbit, the principal tidal frequency X2200 becomes 
nil, for which reason insertion of 

Q ~ , « > (25) 
into (12^ seems to entail an infinitely large torque at the instant of crossing: 

r~ — ~ — ^00, for Y— i>0 , (26) 

a clearly unphysical result. 

Another, very similar objection to ( l23i) originates from the fact that the quality factor is 
inversely proportional to the phase shift: Q ~ 1/e. As the shift ( IT2|) vanishes on crossing 
the synchronous orbit, one may think that the value of the quality factor must, effectively, 
approach infinity. On the other hand, the principal tidal frequency vanishes on crossing the 
synchronous orbit, for which reason (I23p makes the quality factor vanish. Thus we come to a 
contradiction. 

For these reasons, the long-entrenched opinion is that these models introduce discontinuities 
into the expression for the torque, and can thus be considered as unrealistic. 

It is indeed true that, while law fl23|) works over scales shorter than the Maxwell time (about 
10^ yr for most minerals), it remains subject to discussion in regard to longer timescales. 
Nonetheless, it should be clearly emphasised that the infinities emerging at the synchronous- 
orbit crossing can in no way disprove any kind of rheological model. They can only disprove 
the flawed mathematics whence they provene. 

5.2 A case for reasonable doubt 

To evaluate the physical merit of the alleged infinite-torque "paradox" , recall the definition of 
the quality factor. As part and parcel of the linearity approximation, the overall damping inside 
a body is expanded in a sum of attenuation rates corresponding to each periodic disturbance: 

{E) = Y.^^{x.d) (27) 



For the dependence of Af;^ upon Ximpq Efroimsky & Lainey 2007. 
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where, at each frequency Xi j 



( E(y ) ) E iy.) 

( . . . ) designating an average over a flexure cycle, E{xJ denoting the energy of deformation 
at the frequency Xi ? QiXi) being the quality factor of the medium at this frequency. 

This definition by itself leaves enough room for doubt in the above "paradox" . As can be 
seen from ( 128|) . the dissipation rate is proportional not to 1/Q{x) but to x/Qix) ■ This way, 
for the dependence Q ~ 5 the dissipation rate {E) will behave as • the limit of 

X ^ , this scaling law portends no visible difficulties, at least for the values of a up to unity. 
While raising a above unity may indeed be problematic, there seem to be no fundamental 
obstacle to having materials with positive a taking values up to unity. So far, such values of 
a have caused no paradoxes, and there seems to be no reason for any infinities to show up. 



5.3 The phase shift and the quaUty factor 

As another preparatory step, we recall that, rigorously speaking, the torque is proportional 
not to the phase shift e itself but to sin e . From fl?I]) and (1251) we obtain: 

11 

sine = , = = . 29 

i + Q + 8 X 

We see that only for large values of Q one can approximate | sine | with 1/Q (crossing of 
the synchronous orbit not being the case). Generally, in any expression for the torque, the 
factor 1/Q must always be replaced with 1/ a/1 + • Thus instead of f l24l) we must write: 



11 
sme = , = , = , (30) 

/, ^2 I 2a 2a 

1 + Q a/i + ^ X 



£ being a dimensional constant from ( 123|) . 

Though this immediately spares us from the fake infinities at x , we still are facing 
this strange situation: it follows from (!29|) that, for a positive a and vanishing x 1 Ihe phase 
lag e must be approaching 7r/2, thereby infiating the torque to its maximal value (while on 
physical grounds the torque should vanish for zero x)- Evidently, some important details are 
still missing from the picture. 



5.4 The stone rejected by the builders 

To find the missing link, recall that Kaula (1964) described tidal damping by employing the 
method suggested by Darwin (1880): he accounted for attenuation by merely adding a phase 
shift to every harmonic involved - an empirical approach intended to make up for the lack of 
a consistent hydrodynamical treatment with viscosity included. It should be said, however, 
that prior to the work of 1880 Darwin had published a less known article (Darwin 1879), in 
which he attempted to construct a self-consistent theory, one based on the viscosity factor of 
the mantle, and not on empirical phase shifts inserted by hand. Darwin's conclusions of 1879 
were summarised and explained in a more general mathematical setting by Alexander (1973). 
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The pivotal result of the self-consistent hydrodynamical study is the following. When a 
variation of the potential of a tidally disturbed planet, U{r) , is expanded over the Legendre 
functions Pim(sin0) , each term of this expansion will acquire not only a phase lag but also 
a factor describing a change in amplitude. This forgotten factor, derived by Darwin (1879), 
is nothing else but cose . Its emergence should in no way be surprising if we recall that the 
damped, forced harmonic oscillator 

X + 2jx + ulx = F e'^* (31) 

evolves as 

^ {-7 + iV'^?-A2 )t (-7 + iV'^?-A2 ) t F cose i{Xt-e) 

x(t) = Ci e + C2 e + — e , (32) 

where the phase lag is 

tane = 27 A (cj^ - A^) , (33) 

and the first two terms in fl32l) are damped away in time0 

In the works by Darwin's successors, the allegedly irrelevant factor of cos e fell through the 
cracks, because the lag was always asserted to be small. In reality, though, each term in the 
Fourier expansions (ITTl) . ([I6]-[T9]), and (l22l) should be amended with cose^^^^ . For the same 
reason, instead of (l30ll . we should write down: 



a a 



cos e sm e I = , / = , (34) 



1 + Vi + 1 + ^ X 

At this point, it would be tempting to conclude that, since (71) vanishes in the limit of x — ^ , 
for any sign of a , then no paradoxes happens on the satellite's crossing the synchronous orbit. 
Sadly, this straightforward logic would be too simplistic. 

In fact, prior to saying that cose sine — 0, we must take into consideration one more 
subtlety missed so far. As demonstrated in the Appendix, taking the limit of g — > is a 
nontrivial procedure, because at small values of Q the interconnection between the lag and 
the Q factor becomes very different from the conventional Q = cot |e|. A laborious calculation 
shows that, for Q < 1 — it/A , the relation becomes: 



sine cose = ± (3g) 



1/3 



l_4^3Q)2/3 + 0(g4/3) 

5 



which indeed vanishes for Q ^ 0. Both e22oo the appropriate component of the torque 
change their sign on the satellite crossing the synchronous orbit. 

So the main conclusion remains in force: nothing wrong happens on crossing the syn- 
chronous orbit, Q.E.D. 



As demonstrated by Alexander (1973), this example indeed has relevance to the hydrodynamical theory 
of Darwin, and is not a mere illustration. Alexander (1973) also explained that the emergence of the cose 
factor is generic. (Darwin (1879) had obtained it in the simple case of Z = 2 and for a special value of the 
Love number: ki = 1.5 .) 

A further investigation of this issue was undertaken in a comprehensive work by Churkin (1998), which 
unfortunately has never been published in English because of a tragic death of its Author. In this preprint, 
Churkin explored the frequency-dependence of both the Love number and the quality factor within a broad 
variety of rheological models, including those of Maxwell and Voight. It follows from Churkin's formulae that 
within the Voight model the dynamical ^2 relates to the static one as cos e . In the Maxwell and other models, 
the ratio approaches cose in the low-frequency limit. 
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6 Conclusions 



In the article thus far we have punctihously spelled out some assumptions that often re- 
main implicit, and brought to light those steps in calculations, which are often omitted as 
"self-evident" . This has helped us to explain that no "paradoxes" ensue from the frequency- 
dependence ! a = 0.3±0.1, which is in fact the actual dependence found for the 
mantle and crust. 

This preprint is a pilot paper. A more comprehensive treatise on tidal torques is to be 
published. (Efroimsky & Williams 2009) 
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Appendix. 

The lag and the quaUty factor: is the formula Q = cot |e| universal? 

The interrelation between the quality factor Q and the phase lag e is long-known to be 

Q = cot |e| , (35) 

and its derivation can be found in many books. In Appendix A2 of Efroimsky & Lainey(2007), 
that derivation is reproduced, with several details that are normally omitted in the literature. 
Among other things, we pointed out that the interrelation has exactly the form fl35l) only in 
the limit of small lags. For large phase lags, the form of this relation will change considerably. 

Since in section 9 of the current paper we address the case of large lags, it would be worth 
reconsidering the derivation presented in Efroimsky & Lainey (2007), and correcting a subtle 
omission made there. Before writing formulae, let us recall that, at each frequency x the 
spectrum of the deformation, the quality factor (divided by 2 tt ) is defined as the peak energy 
stored in the system divided by the energy damped over a cycle of flexure: 

Qix) - - ^"^"ff , (36) 

where AEcydeix) < as we are talking about energy losses0 

We are considering flexure in the linear approximation. Thus at each frequency x the appropriate energy 
loss over a cycle, AEcycieix) i depends solely on the maximal energy stored at that same frequency, Epgakix) ■ 
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An attempt to consider large lags (all the way up to |e| = vr/2) sets the values of Q/2'k 
below unity. As the dissipated energy cannot exceed the energy stored in a free oscillator, the 
question becomes whether the values of Q/2tt can be that small. To understand that they 
can, recall that in this situation we are considering an oscillator, which is not free but is driven 
(and is overdamped). The quality factor being much less than unity simply implies that the 
eigenfrequencies get damped away during less than one oscillation. Nonetheless, motion goes 
on due to the driving force. 

Now let us switch to the specific context of tides. To begin with, let us recall that the 
dissipation rate in a tidally distorted primary is well approximated by the work that the 
secondary carries out to deform the primary: 

E = - J p V ■ VW d^x (37) 

p , V , and W denoting the density, velocity, and tidal potential inside the primary. The 
expression on the right-hand side can be transformed by means of the formula 

p V ■ VW = V-(pVW)-WV-Vp-WV-{pV)=V-ipVW)-WV-Vp + W^ 

at 

where the W V- Vp and dp/dt terms may be omitted under the assumption that the primary 
is homogeneous and incompressible. In this approximation, the attenuation rate becomes 
simply 



E = - J W ■ (p V W)d-'x = - J pW V ■ n dA , (39) 

n being the outward normal to the surface of the primary, and dA being an element of 
the surface area. It is now clear that, under the said assertions, it is sufficient to take into 
account only the radial elevation rate, not the horizontal distortion. This way, formula ( 137|) . 
in application to a unit mass, will get simplified to 

C standing for the vertical displacement (which is, of course, delayed in time, compared to 
W ) . The amount of energy dissipated over a time interval {to , t) is then 

We shall consider the simple case of an equatorial moon on a circular orbit. At each point 
of the planet, the variable part of the tidal potential produced by this moon will read 

W = Wo cosxt , (42) 

the tidal frequency being given by 

X = 2\n - up\ . (43) 
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Let g denote the surface free-fall acceleration. An element of the planet's surface lying beneath 
the satellite's trajectory will then experience a vertical elevation of 

C = h^— cos(xt - |e|) , (44) 



h2 being the corresponding Love number, and |e| being the posztzf^lj phase lag, which for the 
principal tidal frequency is simply the double geometric angle 5 subtended at the primary's 
centre between the directions to the secondary and to the main bulge: 

\e\ = 2 5 . (45) 

Accordingly, the vertical velocity of this element of the planet's surface will amount to 

u = C, = — h2 X — ~ sin(xt — |e|) = — /i2 X — ~ (sinx^ cos|e| — cosx^ sin|e|) . (46) 

The expression for the velocity has such a simple form because in this case the instantaneous 
frequency x is constant. The satellite generates two bulges - on the facing and opposite sides 
of the planet - so each point of the surface is uplifted twice through a cycle. This entails 
the factor of two in the expression ( H3|) for the frequency. The phase in ( l45l) . too, is doubled, 
though the necessity of this is less evident, - see footnote 4 in Appendix Al to Efroimsky & 
Lainey (2007). 

The energy dissipated over a time cycle T = 27t/x, per unit mass, will, in neglect of 
horizontal displacements, be 

r f dW\ , / , Wo\ dWo f'=^ , , , , , 

l\E , = u [ — ] at = — — n.2 X — — — / cos yt (sm vt cos e — cos xt sm e ) at 

Jo \ dr J V g / A=o ^ ^ ^ II ^ 11^ 

= — fi2 X — sm |e| — / cos xt a[xt) = — /12 7: — tt sm |e| , 



g dr X Jxt = o g 

while the peak energy stored in the system during the cycle will read: 

^peak — I ( ~ ] at = ~ [ — h2 X —~ I ~ — / cosx^ (smxt cos |e| — cosx^ sm |e|) at 



Alx \ dr J V g J dr jt=|.|/x 



cosx^ smxt d[xt) / cos xt d[xt) 

--\e\ X Jxt=\e\ 

In the appropriate expression in Appendix Al to Efroimsky & Lainey (2007), the lower limit of 
integration was erroneously set to be zero. To understand that in reality integration over xt 

12 Were we not considering the simple case of a circular orbit, then, rigorously speaking, the expression for 
W would read not as Wo cosxt but as Wo cosw^j^^jf, the tidal frequency uj^-^^^ taking both positive and 
negative values, and the physical frequency of flexure being x = I'^udai I ■ Accordingly, the expression for ^ 
would contain not cos{xt — |e|) but cos{uj^^^^it — e) . As we saw in equation (24), the sign of e is always the 
same as that of uj^i^^i ■ For this reason, one may simply deal with the physical frequency x = l^udai I 
with the absolute value of the phase lag, |e| . 



, Wo dWo cos lei 

Xh2 ^ 

g or X 
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should begin from |e| , one should superimpose the plots of the two functions involved, cosxt 
and sm{xt — |e|) . The maximal energy gets stored in the system after integration through 
the entire interval over which both functions have the same sign. Hence = |e| as the lower 
limit. 

Evaluation of the integrals entails: 



'peak 



whence 



h,— 



Wo dWo 



g dr 



- COS|6| 



1 / TT 

2 V 2 



sm e 



cycle 



2txE 



peak 





TT 


sin e 




2 cos e 


1 

~ 2 


(f-H' 


j sin e 



tan 


e 






e 


) 


tan 


e 



As can be seen from ( 1501) . both the product sine cose and the appropriate component of the 
torque attain their maxima when Q = 1 — tt / A . 

Usually, |e| is small, and we arrive at the customary expression 



= tan|e| + 0(e^ 



(51) 



In the opposite situation, when Q ^ and |e| 7r/2 , it is convenient to consider the small 
difference 



c - ^ II 

in terms whereof the inverse quality factor will read: 
cot ^ 1 1 



1—1^ cot^ tan,^ — ^ z — arctan^; 



3 ^ 



z^ + 0(z^ 



(52) 



(53) 



where z = tan^ and, accordingly, ^ = arctanz = z —^z^ + -^z^ + 0{z'^) . Formula 
( l53l) may, of course, be rewritten as 



1 - ^z' + 0{z') 
5 



3Q 



(54) 



or, the same, as 



(sg) 



1/3 



1 + - z^ + Oiz^) 
5 



(55) 



While the zeroth approximation is simply z = {3QY^^ + 0{Q) , the first iteration gives: 



1 + - (3g)2/3 + o(g^/3) 

5 



tan^ = z = (3g)^/^ 
with q = (3g)^'^^ playing the role of a small parameter. 



1 + 7 + 0{q') 
5 



(56) 
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We now see that the customary relation (15T!) should be substituted, for large lags, i.e., for 
small^ values of Q , with: 



tan lei 



(sg) 



-1/3 



(57) 



The formula for the tidal torque contains a multiplier sin e cos e , whose absolute value can, 
for our purposes, be written down as 



sin |e| cos |e| = cos.^ sin,^ 



whence 



tan,^ 



l + lq^ + 0{q' 



l + tan^^ l + g2[i + 0(g2)] 



(3Q) 



1/3 



sme cose 



±(3Q) 



1/3 



l_^(3g)2/3 + 0(g4/3) 

5 



1 _ ^(3g)2/3 + 0(g4/3) 




(59) 



an expression vanishing for Q ^ . Be mindful that both egggo ^'^d the appropriate compo- 
nent of the torque change their sign on the satellite crossing the synchronous orbit. 
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